METHOD AND APPARATUS FOR BRAIN WAVE FLUCTUATIONS ANALYSIS 



FIELD OF THE INVENTION 

[0001] The present invention relates to diagnosis-oriented medical apparatus, 
and more particularly, to a method for brain wave fluctuations analysis and a relevant 
apparatus thereof. 

BACKGROUND OF THE INVENTION 

[0002] Theoretically, the information of the physiologic and pathologic situation 
and the functional activities of both human's and animal's brains could be obtained 
through a test on the brain wave signals in the same way as electrocardiograms. Today, 
there are two types of signals used in electroencephalographs, analog signals and digital 
signals (DS). However, because the brain wave signals are very weak and complicated, 
even the digital electroencephalographs with higher anti-jamming properties could not 
meet the demand of clinical medical treatment. Therefore, the significance of 
electroencephalogram is far lower than that of electrocardiogram. 

[0003] It is a research achievement of China's aeromedicine to detect and 
analyze the signals of brain wave fluctuations so as to make a judgment in the functioning 
situation and the diseases of the brain, which means much more than what common 
electroencephalogram does in clinical medicine. As a result, it is necessary for the 
medicinal field to get some methods and relevant apparatus to analyze the signals of brain 
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wave fluctuations thoroughly, comprehensively and precisely. The Chinese patent 
ZL96244175.9, titled "encephalofluctuo gram technology (ET)", is designed to collect, 
magnify and regularly sample brain wave signals, digital-analog convert and pre-treat data 
which is then transmitted to PCs. This invention has improved the brain wave signals 
collecting technique to a brand new level. But this apparatus could not analyze the 
collected brain wave signals. Doctors could only get the flexuosity of the brain wave, i.e. 
the electroencephalogram. And how much information useful for diagnosis doctors could 
get from the electroencephalogram depends on the doctors' discrimination in the 
electroencephalogram. Therefore, this technology could still not meet the demand of the 
clinical medicine. 

SUMMARY OF THE INVENTION 

[0004] An object of the present invention is to provide a method to analyze brain 
wave signals with PC technique and a relevant apparatus, which can get a series of 
parameters shown in multiple ways for providing reference in brain functions inspection 
and disease diagnosis. 

[0005] Hereinafter, definitions and explanation of the technical and scientific 
terms and nomenclatures related to the present invention are provided by reference to 
"New Technique on ET-brain Function Research" (published by National Defense 
Industrial Press, 1994), written by Lei MEL 

[0006] S pedigree: the frequency range for brain wave inspection in the present 
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invention is within supra-slow wave with mHz (milli-HZ) as the unit, which is named as 
supra-slow pedigree and S pedigree for short. According to its frequency S pedigree can 
be separately named as SI, S2, S3 etc. For example, the frequency of lmHz corresponds 
to SI, the frequency of 2mHz corresponds to S2, and likewise other spectrum lines of S 
pedigree are named. 

[0007] Fundamental frequency, the several commonly used frequencies such as 1, 
2, 3, 4, 5, 6, 7, 1 land 13mHz which have close relation with neurotransmitter in the brain. 

[0008] Fundamental pedigree: it is the pedigree corresponding to the 
fundamental frequency. 

[0009] Dominant Frequency, dominant spectrum line, optimal value: rank the 
power values of each frequency in the corresponding lead from high to low, and the 
frequencies of the anterior n values Dl-Dn are called dominant frequencies. Their 
corresponding spectrum lines are called dominant spectrum lines and the power value of 
the dominant frequency is called optimal value. 

[0010] Optimal frequency, optimal spectrum line: the maximum Dl in the 
dominant frequency of each lead is named as optimal frequency and its corresponding 
spectrum line is optimal spectrum line. 

[0011] Continuum frequency: sometimes, the values of the spectrum lines in the 
dominant frequency zone are continuous, such as 2, 3 and 4mHz, which are called 
continuum frequency 

[0012] Different frequency: spectrum lines from non-resonance frequencies in S 
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pedigree (the frequency value is a prime number larger than 13) and their frequency 
multiplication are called different frequencies. 

[0013] Special frequency: a series of frequencies such as 23 mHz, 27mHz, 
28mHz and 29mHz and their harmonic frequencies such as 46rnHz and 54mHz and so on 
in S pedigree are special frequencies. 

[0014] Characteristic lines: lines corresponding to continuum frequency, 
different frequency, special frequency and optimal frequency are generally called 
characteristic lines. 

[0015] A/P, L/R: calculating the ratio of the anterior lead power to the posterior 
power of each frequency according to the lead space distribution positions such as lead F3 
to lead C3, which is called as anterior-to-posterior ratio and written as A/P. Comparing the 
power values of the leads which are at the same position in the right brain and left brain 
and calculating the ratio thereof such as lead F3 to lead F4, which is called as left-to-right 
ratio and written as L/R,. 

[0016] Fluctuation value: generally refers to the values of each pedigree obtained 
during the analysis of S pedigree. 

[0017] The present invention is realized according to the following technical 
solution. 

[0018] In this invention, the brain wave signals are first divided into subsections 
in a certain time interval with computer. Power spectrum of each subsection is then 
analyzed by selecting the maximal power amplitude within 0.5 and 50 Hz, and then 
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several times of power spectrum analyses and frequency spectrum analyses are performed, 
and then power spectrum fluctuations within the supra-slow wave can be acquired, after 
that a series of analyses are performed on said fluctuations and a series of data 
parameters can be acquired which are shown in values, graphs and curves finally. 

[0019] In this invention, the analysis of the fluctuations of brain wave includes at 
least the conventional analysis of brain wave. Also, it could include the analysis of 
fluctuation signals of the power of brain wave, analysis of brain wave fluctuations and S 
pedigree which are performed in turn. Error treatment can also be used to correct the 
spectrum. 

[0020] In this invention, the analysis method of fluctuation of brain wave signals 
involves the sampling of brain wave signal as well. Any lead system or lead combination 
can by used for sampling and electrodes arranging, and the 12 leads from the international 
standard lead system is preferred, with the optimal electrodes positions of F3, F4, C3, C4, 
P3, P4, Ol, 02, F7, F8, T5, T6. 

[0021] The detailed procedures and arithmetic for each analysis are described 
respectively as the following. 

[0022] The said conventional power spectrum analysis includes the following 
procedures. 

[0023] ( 1 ) Analyze the power spectrum of the brain wave of a given time 

length., i.e. to carry out Fourier Transformation directly to the observation data having N 

sampling points x ^^of brain wave x ^ , with the result of X ^ \ 
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[0024] (2) Make a square of its amplitude which is then divided by N. As the 
estimate value of the real power spectrum p ( eJa> ) of x ^ , the power spectrum estimated 

P{e») = ±\X N {a>f 

with cyclogram can be expressed as 

[0025] (3) Calculating X n(°>) by using the fast Fourier Transformation, 

N-l -jllL 

which includes X N (k) = J^x(n)W N nk ,k = 0X...,N -l,W M =e N . And the power 

«=0 

P(k) = ±-\X„(ktf 

spectrum can be expressed as 

[0026] The said analysis of the brain wave power fluctuation signals includes the 

following procedures: 

[0027] ( 1 ) Select window function with small side lobe amplitude and fast 

2tzw 

attenuation, which is expressed as co(ri) = 0.5 -0.5cos(-^-),n = 0,l,..JV-l.to cut the 
brain wave signal x(n) and get x N (n) = x(n)a>(n) . 

[0028] (2) Analyze the power spectrum of the above signal x N (n) = x{n)o(n) 
with the sampling time of T seconds. By FFT calculation, the frequency-domain resolution 
of the frequency spectrum is 1/T. Then, select the power amplitude P with the maximal 
energy and its corresponding frequency f. 

[0029] ( 3 ) Divide the time domain brain wave signal with the total time length 
of N seconds into subsections with the interval of T seconds in time sequence. Then carry 
out the said conventional power spectrum analysis and fluctuation signal analysis of brain 
wave power in order to get the fluctuation signal p(n) with the maximal power amplitude 
and its corresponding frequency fluctuation signals f(n),n=l n, n=N/T. This shows the 
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fluctuation process of the power and its corresponding maximal amplitude during the 
whole period. 

[0030] The said analysis of the fluctuations of the brain wave includes the 
following procedures: 

[0031] ( 1 ) Analyze the fluctuation signals p(n) of the maximal power 
amplitude with the length of n points. 

[0032] ( 2 ) Multiply it with Harming window with the length of n and then carry 
out power spectrum analysis. Since the unit of the time length is N second, frequency 
domain resolution is thus 1/N Hz. Select the spectrum lines within a certain frequencies in 
the analysis result of the power spectrum to form the fluctuations of brain wave. 

[0033] ( 3 ) If the whole sampling time is more than N second, divide it into 
subsection with N seconds as the unit and then carry out the said conventional power 
spectrum analysis, analysis of the fluctuation signals of brain wave power and fluctuations 
of the brain wave for each subsection. 

[0034] The said S pedigree analysis includes the following procedures: 

[0035] (1) Find out several dominant lines Dl-Dn with the maximal amplitude 
from the fluctuations of brain wave of each lead and arrange them by a descending 
order , thus S pedigree of single lead data can be obtained which contains n values. 

[0036] (2) For all N leads, there are N*n dominant spectrum lines. Add 
together all the dominant spectrum lines with the same frequency so as to get the general 
pedigree of the S pedigree. 
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[0037] The above analysis method is called as the first level analysis method in 
the series methods of fluctuation signal analysis of brain wave in this invention. Based on 
the first level analysis method , several further analyses which are called the second level 
analysis can be carried out, including the following 24 items. For the convenience of 
description, each item is marked with a serial number. 

[0038] The analyses based on the analysis of fluctuation signals of brain wave 
power include: (l)calculation of entropy (2) single-frequency competition analysis 

[0039] The analyses based on the S pedigree analysis include: (3)analysis of the 
general pedigree of S pedigree, (4)analysis of fundamental pedigree, (5)analysis of 
optical value, (6) analysis of reversion of A/P and analysis of L/R unbalance, (7)analysis 
of special frequency, (8)analysis of different frequency, (9)analysis of continuum frequency, 
(10)analysis of optimal frequency, (11) analysis of space distribution of S pedigree power, 
Q2)analysis of single-frequency power and distribution of its correspondinrelative value 
(L/R) , (13)analysis of the distribution of mean power, (14)analysis of the relative value A/P 
and L/R of power, (15)analysis of the curve of long temporal S pedigree, (16)analysis of the 
curve of long temporal dominant spectrum lines, (©analysis of the curve of long temporal 
fundamental pedigree, (18)analysis of the curve of long temporal space distribution of 
power, (19)analysis of the curve of long temporal entropy, (20)analysis of the curve of long 
temporal special frequency, (21) analysis of the curve of long temporal continuum 
frequency, (22) analysis of the curve of long temporal space distribution of fundamental 
pedigree, (23) analysis of the curve of long temporal conventional power spectrum, (24) 
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distinguishing of long temporal event mark. 

[0040] In the above four items of (7)analysis of special frequency, (8)analysis of 
different frequency, (9)analysis of continuum frequency and (©analysis of optimal 
frequency can be named as characteristic spectrum analysis. 

[0041] The above three items of (12)analysis of single-frequency power and 
distribution of its relative value (L/R) , (13)analysis of the distribution of mean power and 
(14)analysis of the relative value A/P and L/R of power can be named as analysis of power 
distribution. 

[0042] The above ten items of (15)analysis of the curve of long temporal S 
pedigree, (©analysis of the curve of long temporal dominant spectrum lines, (17)analysis 
of the curve of long temporal fundamental pedigree, (©analysis of the curve of long 
temporal space distribution of power, (19)analysis of the curve of long temporal entropy, 
(©analysis of the curve of long temporal special frequency, (21 ) analysis of the curve of 
long temporal continuum frequency, (22) analysis of the curve of long temporal space 
distribution of fundamental pedigree, ( 23 ) analysis of the curve of long temporal 
conventional power spectrum and (24) distinguishing of long temporal event mark can be 
named as analysis of long temporal dynamic curve. 

[0043] All the analysis methods of second level will be described as following. 

[0044] The said (1) calculation of entropy is carried out based on the analysis of 
the fluctuation signals of brain wave and the method is includes: 

13 

[0045] ( 1 ) Calculate entropy according to H = -^p k lg 2 P k > * n which p k 

k=8 
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is the probability of each frequency being optimal frequency in the brain wave. 

[0046] (2 ) Calculate the total entropy from the probability distribution of all the 
N leads (total is n*N) . 

[0047] The said (2) analysis of the single- frequency competition is also based on 
the analysis of the fluctuation signals of brain wave and the detailed method is to 
accumulate the number of the same optimal frequencies in frequency fluctuations of brain 
wave fluctuation signals along with time changing and so as to get the probability curve of 
the optimal frequency 

[0048] The following analysis methods from No. (3) to ( 24 ) are all based on the 
analysis of S pedigree. 

[0049] The analysis method of the said (3) analysis of the general pedigree of S 
pedigree is to present the data of the general pedigree of S pedigree acquired during S 
pedigree analysis with graphs. 

[0050] The analysis method of the said (4) analysis of the fundamental pedigree 
is to analyze statistically of the S pedigree corresponding to the fundamental frequency in 
S pedigree. Add the values of the multiple periodic frequency thereof from 3 
milli-Hz(mHz) (for example, when 3mHz is analyzed the value of 6mHz, 9mHz and 
etc. should also be added). Meanwhile, all the leads are classified into multiple 
subsections such as front, rear, left, and right subsections according to the positions of the 
leads in the head for respective statistic analysis. 

[0051] The analysis method of the said (5) analysis of the optimal value is to 
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present the power value of the dominant frequency and its corresponding frequency 
according to the space position distribution of the lead. 

[0052] The analysis method of the said (6) analysis of reversion of A/P and 
analysis of L/R unbalance is to calculate the anterior-to-posterior ratio A/P of the power 
value of each frequency according to the space distribution of the lead and then present the 
frequencies whose A/P values are more than a predetermined limited value. Meanwhile, 
calculate the left-to-right ratio L/R and present the frequencies whose L/R values are more 
than a predetermined limited value. 

[0053] The analysis method of the said four items of (7)analysis of special 
frequency, (8)analysis of different frequency, (9)analysis of continuum frequency, (10) 
analysis of optimal frequency is to present the special frequency, different frequency, 
continuum frequency and optimal frequency of each lead according to the space 
distribution position of the lead. 

[0054] The analysis method of the said (IDanalysis of space distribution of S 
pedigree power is to arrange all the power values of each spectrum line in the brain wave 
fluctuations according to their space lead positions. Open a window in the display 
interface in the form of "sub-interface in interface" to select the spectrum lines. For the 
selected spectrum line, its power values of each lead are shown according to the 
space distribution position of the lead. 

[0055] The analysis method of the said (12) analysis of single- frequency power 
and distribution of its relative value (L/R) is to add together all the power corresponding to 
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the dominant spectrum lines Dl-Dn of each lead so as to get the total power value of 
each lead. Then present the power value of the fundamental frequency and the leflt-to-right 
ratio (L/R) which is more than a predetermined limited value or less than the reciprocal of 
said predetermined limited value according to the space distribution position of the lead. 

[0056] The analysis method of the said (13)analysis of the distribution of mean 
power is to present the mean power value of each lead according to the space 
distribution position of the lead. 

[0057] The analysis method of the said (14)analysis of the relative value A/P and 
L/R of power is to calculate the anterior-to-posterior ratio and the left-to-right ratio of the 
power value according to the space distribution of the lead. 

[0058] The analysis method of the said 05)analysis of the curve of long temporal 
S pedigree is to form a curve by using fluctuation value of each pedigree or each spectrum 
line of each lead or all leads as vertical axis and time as horizontal axis. Open a window 
in the display interface in the form of "sub-interface in interface" for spectrum lines or 
pedigree selection. 

[0059] The analysis method of the said (l©analysis of the curve of long temporal 
dominant spectrum lines is to form a curve by using the frequency of the spectrum lines 
within the dominant spectrum line zone as the vertical axis and time as horizontal axis. 
Open a window in the display interface in the form of "sub-interface in interface" for the 
arrangement selection of dominant spectrum lines (Dl-Dn). 

[0060] The analysis method of the said (17)analysis of the curve of long temporal 
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fundamental pedigree is to form a curve by using the fluctuation values in each time 
interval of the fundamental pedigree of all leads or each lead as vertical axis and time as 
horizontal axis. Open a window in the display interface in the form of "sub-interface in 
interface" for pedigree selection. 

[0061] The analysis method of the said (18)analysis of the curve of long temporal 
space distribution of power is to form a curve by using the power value of each spectrum 
line of each lead as vertical axis and time as horizontal axis. Open a window in the 
display interface in the form of "sub-interface in interface" for spectrum line selection. 

[0062] The analysis method of the said (19)analysis of the curve of long temporal 
entropy is to make out a curve by using entropy of all leads or each lead as vertical axis 
and time as horizontal axis to present the changes of entropy along with time. 

[0063] The analysis method of the said (20)analysis of the curve of long temporal 
special frequency is to make out a curve by using the number of special frequencies 
appeared for each lead or all leads as vertical axis and time as horizontal axis so as to 
observe the changes of special frequency along with time. 

[0064] The analysis method of the said (21) analysis of the curve of long 
temporal continuum frequency is to form a curve by using the number of continuum 
frequencies appeared for each lead or all leads as vertical axis respectively and time as 
horizontal axis so as to observe the changes of continuum frequency along with time. 

[0065] The analysis method of the said (22) analysis of the curve of long 
temporal space distribution of fundamental pedigree is to form curves each of which is 
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presented according to the space position distribution of the lead by using the power value 
of each lead which is acquired from the result of analysis of said single- frequency power 
and the distribution of its relative value (L/R) as vertical axis and time as horizontal axis. 
Open a window in the display interface in the form of "sub-interface in interface" for 
pedigree selection. 

[0066] The analysis method of the said (23) analysis of the curve of long 
temporal conventional power spectrum is to form n dynamic curves by using the power 
values of the multiple frequencies Dl-Dn selected from the conventional power spectrum 
which have maximal amplitudes and which are arranged in a descending order as vertical 
axis and time as horizontal axis. 

[0067] The analysis method of the said (24) distinguishing of long temporal 
event mark is to distinguish the event mark signals recorded in a certain brain wave 
recording box and present them in the playback of brain wave signals and in the 
corresponding position of the time axis of various dynamic curves. 

[0068] The above mentioned second level analysis includes 24 items. By 
selecting and combining said 24 items the solutions can be as follows: 

[0069] ( 1 ) all items are used at the same time; 

[0070] (2) except for item of "distinguishing of event mark", select any 
individual item; 

[0071] (3) Use any arbitrary combination of any item except for item of 
"distinguishing of event mark" 
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[0072] (4) Combine the item of distinguishing of event mark with any one or 
ones of the other nine items of the items of long temporal dynamic curve analysis ; 

[0073] The above mentioned methods can be used to treat the brain wave signals 
sampled by any one lead or any combination of more leads and the operation result of any 
one lead or leads combination in the analysis result can also be selected to output to a 
terminal processor to display, print or store. 

[0074] This invention also reveals a kind of apparatus for the analysis of the 
fluctuation signals of brain wave. It includes electrodes, digital brain wave signal 
amplifier or brain wave recorder, PC, data processor and terminal processor. They are 
linked together in sequence. The electrodes are used to sample the brain wave signals 
which are transmitted to the digital brain wave signal amplifier and/or brain wave recorder 
for data acceptance, amplification, digital/analog conversion, digital filter or/and data 
storage. Data in the digital brain wave signal amplifier and/or brain wave recorder will be 
transmitted up to PC and the data treatment and fluctuation analysis will be performed by 
the data processor connected with the PC . And the analytical result will be transmitted 
to a terminal processor for storage, display or printing. 

[0075] Any lead or combination of multiple leads of any lead connection method 
can be selected for placing the said electrodes . 

[0076] The said data processor includes such modules for the first level analysis 
as module for analysis of conventional power spectrum, module for analysis of the 
fluctuation signals of the brain wave, module for the analysis of the fluctuations of brain 
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wave and module for analysis of S pedigree. They are connected in sequence. Data 
output from the former module is input to the latter module for calculation. The data 
processor can also include any module of the following 24 modules for performing the 
second level analysis. The second level analysis module accepts the result input from the 
first level and perform further analysis. For the convenience of description, the second 
analysis modules are all marked with series numbers. 

[0077] Modules which are connected to the module for analysis of the power 
fluctuation signals of brain wave and accept its data include (Dmodule for analysis of 
entropy calculation and (2)module for analysis of single- frequency competition. 

[0078] Modules which are connected to the module for analysis of supra-slow(S) 
pedigree and accept its data include: (3)module for analysis of the general S pedigree, 
(4)module for analysis of fundamental pedigree, (5)module for analysis of optical value, 
(6) module for analysis of reversion of A/P and analysis of L/R unbalance, (7)module for 
analysis of special frequency, (8)module for analysis of different frequency, (9)module for 
analysis of continuum frequency, (10)module for analysis of optimal frequency, (11) module 
for analysis of power space distribution of S pedigree , (12)module for analysis of power of 
single- frequency and distribution of its relative value (L/R), (13)module for analysis of the 
distribution of mean power, Q4)module for analysis of the power relative value of A/P and 
L/R , (15)module for analysis of the curve of a long temporal S pedigree, (16)module for 
analysis of the curve of a long temporal dominant spectrum lines, (17)module for analysis 
of the curve of a long temporal fundamental pedigree, (18)module for analysis of the 
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curve of a long temporal power space distribution r, (19)module for analysis of the curve of 
a long temporal entropy, (20)module for analysis of the curve of a long temporal special 
frequency, ( 21 ) module for analysis of the curve of a long temporal continuum frequency, 
(22) module for analysis of the curve of a long temporal space distribution of 
fundamental pedigree,(23)module for analysis of the curve of a long temporal 
conventional power spectrum,(24)module for distinguishing of long temporal event mark. 

[0079] All modules for both two levels of analyses are described respectively 
hereinafter . Wherein, modules for the first level analysis , 

[0080] The said module for conventional power spectrum analysis is used to 
analyze the power spectrum of the time domain brain wave with the given time period so 
as to get the power spectrum which can be expressed as P{k) =^\X N (k)\ 2 . 

[0081] The said module for analysis of brain wave power fluctuation signals 

are used to select a window function with small side lobe amplitude and fast attenuation, 

which is expressed as coin) = 0.5 - 0.5cos(-^^-),« = 0,1 V .JV-1. Cut the brain wave 

N 

signal x(ri) and get x N (n) = x(n)co{ri) . 

[0082] Analyze the power spectrum of the above mentioned 
signal x N (n) = x{ri)co(n) with T seconds as the sampling time. Processed by FFT , the 
frequency-domain resolution of the frequency spectrum is acquired i.e. 1/T. The power 
amplitude P with the maximal energy and its corresponding frequency f are selected. The 
time domain brain wave signal are divided with the total time length of N seconds into 
subsections with the interval of T seconds in time sequence. Then carry out the said 
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module for conventional power analysis and module for analysis of power fluctuation 
signal of brain wave in order to get the time fluctuation signal of the maximal power 

amplitude p(n) and its corresponding frequency fluctuation signal f(n),n=l n, n=N/T. 

[0083] The said module for analysis of the fluctuations of the brain wave is 
used to analyze the fluctuation signal of maximal power amplitude with the length of n 
points which is expressed as p(n). Multiply said p(n) with Hanning window with the 
length of n and then analysis of power spectrum is carried out. Since the unit of the time 
length is N seconds, frequency domain resolution is 1/N Hz. The spectrum lines within a 
certain frequency zone are selected from the analysis result of the power spectrum to form 
the fluctuations of brain wave. If the whole sampling time is more than N seconds, it will 
be divided into subsections with N seconds as the unit and said module for analysis of 
conventional power spectrum the module for analysis of the power fluctuation signals of 
brain wave and the module for analysis of fluctuations of the brain wave are performed 
for each subsection. 

[0084] The said module for S pedigree analysis is used to find out several 
dominant spectrum lines Dl-Dn with the maximal amplitude from the fluctuations of 
brain wave of each lead which are arranged in a descending order . In this way, S pedigree 
with single lead data can be obtained, totally n values. For all N leads, there are N*n 
dominant lines. By adding together all the dominant 1 lines with the same frequency the 
general pedigree of S pedigree can be obtained. 

[0085] Modules for the second level analysis are described respectively 
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hereinafter . 

[0086] The data acquired from the result of the analysis performed by the 
module for analysis of brain wave power fluctuations signals for first level analysis are 
transmitted to (1) module for calculation of entropy and (2) module for analysis of 
single- frequency competition to perform the second level analysis, wherein, 

[0087] The said module for entropy calculation is used to calculate entropy 
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according to H = -^p k lg 2 p k , in which p k is the probability of each frequency being 

optimal frequency in the brain wave and to calculate the total entropy from the general 
probability distribution of all the N leads (total is n*N) . 

[0088] The said module for single-frequency competition analysis is used to 
accumulate the number of the same optimal frequencies in frequency fluctuations of brain 
wave fluctuation signals along with time changing (1-n data sections) so as to get the 
probability curve of the optimal frequency. 

[0089] Data from module for S pedigree analysis for the first level analysis is 
transmitted to the modules for the second level analysis. 

[0090] The analysis module of the said (3) module for analysis of the general 
pedigree of S pedigree is to present the data of the general pedigree of S pedigree acquired 
during S pedigree analysis with graphs. 

[0091] The analysis module of the said (4) module for analysis of the 
fundamental pedigree is to analyze statistically the S pedigree corresponding to the 
fundamental frequency in S pedigree. Add the values of the multiple periodic frequency 
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thereof when the frequency is higher than 3 mHz, 3mHz is included (for example, when 
3mHz is analyzed , the value of 6mHz, 9mHz and etc. should also be added). Meanwhile, 
all the leads are classified into multiple subsections such as front, rear, left, and right 
subsections according to the positions of the leads in the head for respective statistical 
analysis . 

[0092] The analysis module of the said (5) module for analysis of the optimal 
value is to present the power value of the dominant frequency and its corresponding 
frequency according to the space position distribution of the lead. 

[0093] The analysis module of the said (6) module for analysis of reversion of 
A/P and for analysis of L/R unbalance is to calculate the anterior-to-posterior ratio A/P 
of the power value of each frequency according to the space distribution of the lead and 
then present the frequencies whose A/P values are more than a predetermined limited 
value. Meanwhile, calculate the left-to-right ratio L/R and present the frequencies whose 
L/R values are more than a predetermined limited value. 

[0094] The analysis module of the said four modules which include (7)module 
for analysis of special frequency, (8)module for analysis of different frequency, (9)module 
for analysis of continuum frequency, (10)module for analysis of optimal frequency are used 
to present the special frequency, different frequency, continuum frequency and optimal 
frequency of each lead according to the space distribution position of the lead 
respectively. 

[0095] The analysis module of the said ODmodule for analysis of power space 
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distribution of S pedigree is to arrange all the power values of each spectrum line in the 
brain wave fluctuations according to their space lead positions. A window is opened in the 
display interface in the form of "sub-interface in interface" to select the spectrum lines. 
For the selected spectrum line, r its power values of each lead are shown according to the 
space distribution position of the lead. 

[0096] The analysis module of the said (12) module for analysis of power of 
single-frequency and distribution of its relative value (L/R) is to add together all the power 
values corresponding to the dominant spectrum lines Dl-Dn of each lead so as get the 
total power value of each lead. Then present the power value of the fundamental 
frequency and the left-to-right ratio (L/R) which is more than a predetermined limited 
value or less than the reciprocal of said predetermined limited value according to the space 
distribution position of the lead. 

[0097] The analysis module of the said (13)module for analysis of the distribution 
of mean power is to present the mean power value of each lead according to the space 
distribution position of the lead. 

[0098] The analysis module of the said (14)module for analysis of the power 
relative value A/P and L/R is to calculate the anterior-to-posterior ratio and the left-to-right 
ratio of the power value according to the space distribution of the lead. 

[0099] The analysis module of the said (15)module for analysis of the curve of 
long temporal S pedigree is to form a curve by using fluctuation value of each pedigree or 
each spectrum line of each lead or all leads as vertical axis and time as horizontal axis. A 
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window is opened in the display interface in the form of "sub-interface in interface" for 
spectrum lines or pedigree selection. 

[0100] The analysis module of the said (16)module for analysis of the curve of 
long temporal dominant spectrum lines is to form a curve by using the frequency of the 
spectrum lines within the dominant spectrum line zone as the vertical axis and time as 
horizontal axis. And a window is opened in the display interface in the form of 
"sub-interface in interface" for the arrangement selection of dominant spectrum lines 
(Dl-Dn).The analysis module of the said (17)module for analysis of the curve of long 
temporal fundamental pedigree is to form a curve by using the fluctuation values in each 
time interval of the fundamental pedigree of all leads or each lead as vertical axis and time 
as horizontal axis. And a window is opened in the display interface in the form of 
"sub-interface in interface" for pedigree selection. The analysis module of the said 
(18)module for analysis of the curve of long temporal space distribution of power is to form 
a curve by using the power value of each spectrum line of each lead as vertical axis and 
time as horizontal axis. And a window is opened in the display interface in the form of 
"sub-interface in interface" for spectrum line selection. 

[0101] The analysis module of the said (19)module for analysis of the curve of 
long temporal entropy is to form a curve by using entropy of all leads or each lead as 
vertical axis and time as horizontal axis to present the changes of entropy along with time.. 

[0102] The analysis module of the said (20)module for analysis of the curve of 
long temporal special frequency is to form a curve by using the number of special 
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frequencies appeared for each lead or all leads as vertical axis and time as horizontal axis 
so as to observe the changes of special frequency along with time. 

[0103] The analysis module of the said (21) module for analysis of the curve of 
long temporal continuum frequency is to form a curve by using the number of continuum 
frequencies appeared for each lead or all leads as vertical axis respectively and time as 
horizontal axis so as to observe the changes of continuum frequency along with time. 

[0104] The analysis module of the said (22 ) module for analysis of the curve of 
long temporal space distribution of fundamental pedigree is to form curves each of which 
is presented according to the space position distribution of the lead by using power value 
of each lead which is acquired from the result of analysis of said single-frequency power 
and the distribution of its relative value (L/R) as vertical axis and time as horizontal axis. 
And a window is opened in the display interface in the form of "sub-interface in interface" 
for pedigree selection. 

[0105] The analysis module of the said (23) module for analysis of the curve of 
long temporal conventional power spectrum is to form n dynamic curves by using the 
power values of the multiple frequencies Dl-Dn selected from the conventional power 
spectrum which have maximal amplitude and which are arranged in a descending order as 
vertical axis and time as horizontal axis. 

[0106] The analysis module of the said (24) module for distinguishing of long 
temporal event mark is used to recognize the event mark signals recorded in a certain 
brain wave recorder and present them in the playback of brain wave signals and present 
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them in the corresponding position of the time axis of various dynamic curves. It can be 
connected to each long temporal analysis module respectively. 

[0107] Data of brain wave signal of any lead or any combination of multiple 
leads sampled by the electrodes can be selected for treating by the said data processor. 

[0108] The above mentioned modules for the second level analysis can be 
performed at the same time or it is also applicable to perform any combination of the said 
modules selected from said module (1) ~ module (23) except for module (24) i.e. module 
for distinguishing of event mark. 

[0109] The above-mentioned analysis methods can be used to treat data sampled 
by one lead or combination of multiple leads. The operation result of one lead or 
combination of multiple leads can also be selected from the analysis result and be output 
to the terminal processor. 

[0110] The said terminal processor includes display, printer and storage 
apparatus such as hardware disk, floppy disk and compact disk etc, which is used to accept 
the signals from data processor and select the operation result of one lead or combination 
of more leads from data processor to store, display and print. 

[0111] According to the method and apparatus of the present invention, the data 
and curve displayed which has been inspected and calculated can be used to analyze the 
situation of the brain function and neurotransmitter in the patients' brain, and can also be 
used to analyze the functional changes of the brain which is impossible for CT and 
NMR. It provides direct and objective index for the diagnosis of functional 
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encephalopathy in the medical field and make up the blank of objective inspective index 
in functional encephalopathy such as lunacy in the medical field. 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0112] Fig. 1 is a block diagram according to a first embodiment of the present 
invention; 

[0113] Fig. 2 is a block diagram according to a second embodiment of the 
present invention; 

[0114] Fig. 3 is a block diagram according to a third embodiment of the present 
invention; 

[0115] Fig. 4 is a block diagram showing structures of a data processor of the 
present invention. 

[0116] Wherein, if the blocks are connected by arrowhead, it means logic 
relations. In this condition, data acquired in the anterior module is the basis of the 
operation of the latter module. That's to say, the latter module relies on the result of the 
anterior module. If the blocks are connected by beeline, it means inclusion relation. In this 
case, modules in the anterior blocks are composed of the latter modules. If the anterior 
block are connected by dashed, it means the anterior shows only the general name of the 
sort of the latter modules' attributes. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
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[0117] Further description is given in the combination of the drawings 
hereinafter. 

[0118] The process of analyzing the brain wave fluctuations includes: getting 
the brain wave fluctuation signals from the brain wave with computer fluctuation scanning 
technique, then analyzing of the frequency spectrum of the fluctuation signals in order to 
get the power spectrum within the scope of 1 to 255 mHz. The detailed analysis process 
can be summarized as: dividing the 1024 seconds brain wave data into 512 subsections 
each of which is with 2 seconds , analyzing the power spectrum of data of each subsection 
and select the maximal power amplitude within the scope of 0.5 to 50Hz, analyzing the 
power spectrum of the time fluctuations of the maximal power amplitude (that is the 
maximal values of the corresponding power spectrums of said 521 subsections) so as to 
get the fluctuations of the power spectrum within the scope of 1 to 255mHz. 

[0119] There are three types of structure of the apparatus according to the 
present invention, which are as follows. 

[0120] As shown in Fig. 1, it includes electrodes A, a digitized brain wave 
amplifier B, a USB interface J, a Personal Computer C, a data processor D and a terminal 
processor E, which are connected in turn. 

[0121] As shown in Fig. 2, it includes electrodes A, a brain wave recorder F, a 
USB interface J, a PC C, a data processor D and a terminal processor E, which are 
connected in turn. 

[0122] As shown in Fig. 3, it includes an electrode A, a digital brain wave 
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amplifier B and a brain wave recorder F, a USB interface J, a PC C, a data processor D 
and a terminal processor E, which are connected in turn. Wherein, the electrode A is 
connected to the digital brain wave amplifier B and the brain wave recorder at the same 
time. The digital brain wave signal amplifier B is connected to the PC C, and the brain 
wave recorder F is connected to the PC C through the USB interface J. Finally, the PC C is 
connected to the data processor D and the terminal processor E in sequence. 

[0123] The functions of each component in the present invention are described 
below. 

[0124] The electrode A is used to sample brain wave signals. Location of the 
electrode complies with the 12 leads international standard lead system with the position 
separately in F3^ F4^ C3> C4> P3> P4s OK F7> F8^ T5> T6. Frequency of brain 
sampling is 128Hz. 

[0125] The digital brain wave signal amplifier B has functions of accepting 
signals, signal amplification, digital/analog conversion and data filter etc. 

[0126] The brain wave recorder F is used to sample and analyze long temporal 
(longer than 18min) brain wave data. It has functions of signal sampling, signal 
amplification, digital/analog conversion and data filter, data storage and data playback etc. 
Data in the brain wave recorder is uploaded to PC for fluctuation analysis. 

[0127] The PC C is configured to have a PIV host and 256M EMS memory. 

[0128] The terminal processor E is composed of display, printer and storage 
devices such as hardware disk, floppy disk and compact disk. It accepts data from the 
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data processor D for display, storage or print. 

[0129] The data processor D includes an analysis module of module for a first 
level conventional power spectrum power Dl, an analysis module of fluctuation signals of 
brain wave D2, an analysis module of fluctuation of brain wave D3 and an analysis 
module of S pedigree D4. They are connected to each other in sequence. And data 
output from the anterior module is transmitted to the latter module for analysis. The data 
processor can also include any module of the 24 modules for the second level analysis. 
For the convenience of description, the second level analysis modules are all marked with 
series number. 

[0130] Modules which are connected to module D2 and accept its data are : 
(1) an analysis module of entropy calculation D5 and (2) an analysis module of 
single- frequency competition D6. 

[0131] Modules which are connected to module D4 and accept its data are : 
(3)an analysis module of the general pedigree of S pedigree (i.e. module for analysis of the 
general pedigree of S pedigree, and for the other modules )D4a, (4) an analysis module of 
fundamental pedigree D4b, (5)an analysis module of optical value D4c, (6)an analysis 
module of the reversion of A/P; analysis module on L/R unbalance D4d, (7)an analysis 
module of special frequency D4ea, (8)an analysis module of different frequency D4eb, 
(9)an analysis module of continuum frequency D4ec, (10)an analysis module of optimal 
frequency D4ed, (ll)an analysis module of space distribution of S pedigree power D4f, 
(12)an analysis module of power distribution of single- frequency D4ga, (13)an analysis 
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module of the distribution of mean power D4gb, (14)an analysis module of the relative 
value A/P and L/R of power D4gc, (15)an analysis module on the curve of long temporal S 
pedigree D4ha, (16)an analysis module of the curve of long temporal optical lines D4hb, 
Q7)an analysis module of the curve of long temporal fundamental pedigree D4hc, (18)an 
analysis module of the curve of long temporal space distribution of power D4hd, (19)an 
analysis module of the curve of long temporal entropy D4he, (20)an analysis module of 
the curve of long temporal special frequency D4hf, (21) an analysis module of the 
curve of long temporal continuum frequency D4hg, ( 22 ) an analysis module of the curve 
of long temporal space distribution of fundamental pedigree D4hh, ( 23 ) an analysis 
module of the curve of long temporal conventional power spectrum D4hi, (24) a 
distinguishing module of long temporal event mark D4hj. 

[0132] Based on S pedigree analysis, (3)analysis of the general pedigree of S 
pedigree, (4)analysis of fundamental pedigree, (5)analysis of optical value, (6) reversion of 
A/P; analysis of L/R unbalance, (7)analysis of special frequency, (8)analysis of different 
frequency, (9)analysis of continuum frequency, (10)analysis on optimal frequency, (11) 
analysis of space distribution of S pedigree power, (12)power of single-frequency and 
distribution of its corresponding value (L/R), (13)analysis of the distribution of mean power, 
(14)analysis of the relative value A/P and L/R of power, (15)analysis of the curve of long 
temporal S pedigree, (16)analysis of the curve of long temporal optical lines, (17)analysis of 
the curve of long temporal fundamental pedigree, (18)analysis of the curve of long temporal 
space distribution of power, (19)analysis of the curve of long temporal entropy, (20)analysis 

29 



of the curve of long temporal special frequency, (21) analysis of the curve of long 
temporal continuum frequency, (22) analysis of the curve of long temporal space 
distribution of fundamental pedigree, ( 23 ) analysis of the curve of long temporal 
conventional power spectrum, (24) distinguishing of long temporal event mark. 

[0133] Wherein, the above (7)analysis module of special frequency D4ea, 
(8)analysis module of different frequency D4eb> (9)analysis module of continuum 
frequency D4ec, (©analysis module of optimal frequency D4ed are integrated as an 
analysis module of characteristic lines D4e. 

[0134] The above (12) analysis module of power distribution of single-frequency 
D4ga, (13)analysis module of the distribution of mean power D4gb, (14)analysis module of 
the relative value A/P and L/R of power D4gc are integrated as an analysis module of 
power distribution D4g. 

[0135] The above ten modules including Q5)the analysis module of the curve of 
long temporal S pedigree D4ha, Q6)the analysis module of the curve of long temporal 
optical lines D4hb, (17)the analysis module of the curve of long temporal fundamental 
pedigree D4hc, (18)the analysis module of the curve of long temporal space distribution of 
power D4hd, (19)the analysis module of the curve of long temporal entropy D4he, (20)the 
analysis module of the curve of long temporal special frequency D4hf, (21) the analysis 
module of the curve of long temporal continuum frequency D4hg, (22) the analysis 
module of the curve of long temporal space distribution of fundamental pedigree D4hh, 
( 23 ) the analysis module of the curve of long temporal conventional power spectrum 
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D4hi, (24) the distinguishing module of long temporal event mark D4hj are integrated as 
an analysis module of long temporal dynamic curve D4h. 

[0136] The functions of each module are described as below. 

[0137] 1, Analysis module of conventional power spectrum Dl 

[0138] Carry out Fourier Transformation directly to the observation data 
which has N sampling points of brain wave x ^ n \ with the result ofX N (e ja> ). Make a 
square of its amplitude which is then divided by N. The result is the estimate value of the 

real power spectrum p ^ eJ ) of x ^ n ^ . The power spectrum estimated with cyclogram can 
be expressed as P(e ja> ) = — \X N (co)\ 2 . 

[0139] Calculation of X N (co) is acquired by using fast Fourier Transformation 

X N {k) = ^x(ji)W N nk ,k - 0,1,..., N — l,W N =e N . And the power spectrum can be 

77 = 0 

expressed as iV 

[0140] Analyze the power spectrum of a wave signal in the time domain and the 

given time length is 8s so as to get the energy distribution in the frequency domain i.e. 

conventional power spectrum. Then transmit the power spectrum to the terminal 

processor E in the form of graph. 

[0141] 2. Analysis module of fluctuation signals of brain wave power D2 
[0142] During process of the real power estimation practice, the window 

function with small major side lobe amplitude and fast attenuation should be selected. The 

major lobe of Harming window is relatively wide, but it has small side lobe and high 
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attenuation speed. It can be expressed as a>(ri) = 0.5-0.5cos( ),n = 0,1,..JV-1. . 

Cutting the data with Hanning window, the result is x N (n) = x(ri)co{ri) . 

[0143] Multiplying the time domain brain wave signal of a given length of 2s 
with Hanning window, and analyze the power spectrum of the product. . Since the 
sampling time is 2s, perform operation of FFT and the frequency resolution of the 
frequency spectrum is obtained, which is l/T=l/2=0.5Hz (T is the sampling time). Select 
the power amplitude p with maximal energy and its corresponding frequency f within the 
frequency zone from 8Hz to 13Hz (8Hz, 8.5Hz, 9Hz, 9.5Hz, 10Hz, 10.5Hz, 11Hz, 
11.5Hz, 12Hz, 12.5Hz, 13Hz, totally including 11 values ) . 

[0144] Divide the 1024 seconds brain wave data in time sequence into 512 
subsections each of which is of 2 seconds . Perform the analysis module of conventional 
power spectrum Dl and analysis module of the brain wave fluctuation signals D2 to get 
the time fluctuation signals of the maximal power amplitude p(n) and its corresponding 

frequency fluctuation signals f(n), n=l 512. It presents the fluctuation process of the 

power and its corresponding maximal amplitude within 1024 seconds. 

[0145] Transmit the time fluctuation signals to analysis module of fluctuation 
signals of brain wave D3, analysis module of entropy calculation D5 and analysis module 
of single-frequency competition D6. And transmit it to the terminal processor E at the 
same time. 

[0146] 3 . Analysis module of brain wave fluctuations D3 

[0147] Analyze the fluctuation signals p(n) of the maximal power amplitude 
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within the length of 1024s which is of 512 points. Multiply it with Hanning window with 
the length of 5 12 and then carry out the analysis of power spectrum. Since the unit of the 
time length is 1024 seconds, the frequency domain resolution is 1/1024 Hz. The spectrum 
lines within the frequencies of 1/I024*(lto255) Hz from the analysis result of the power 
spectrum form the fluctuations of brain wave. If the whole sampling time is more than 
1024 seconds, divide it into subsection with 1024 seconds as the unit and then carry out 
the said analysis module of conventional power spectrum Dl, analysis module of the 
fluctuation signals of brain wave power D2 and analysis modules of fluctuations of the 
brain wave D3 for each subsection repeatedly. 

[0148] The analysis module of brain wave fluctuations D3 will transmit the 
fluctuations signals of brain wave to the analysis module of S pedigree D4 and also to the 
terminal processor E. 

[0149] 4 . Analysis module of S pedigree D4 

[0150] Find out 8 dominant lines D1-D8 with maximal amplitude from the 
fluctuations of brain wave of each lead and arrange them by a descending order so as to 
get the single lead S pedigree with 8 values. And for all 12 leads, there are 12x8=96 
dominant spectrum lines. Add together the number of dominant spectrum lines with the 
same frequency, and then the general pedigree of S pedigree is obtained. 

[0151] Analysis module of S pedigree D4 transmits the signals of single-lead S 
pedigree and general S pedigree to analysis module on the general peigree of S pedigree 
D4a, analysis module of fundamental pedigree D4b, analysis module of optical value D 
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(1-8) D4c, module for analysis of the reversion of A/P and L/R unbalance D4d, analysis 
module of characteristic lines D4e, analysis module of space distribution of pedigree 
power D4f, analysis module of power distribution D4g and analysis module of long 
temporal dynamic curve D4h. It also transmits signals of general S pedigree to the 
terminal processor E. 

[01 52] 5 . Entropy calculation module D5 

13 

[0153] Calculate entropy according to H = -^p k lg 2 p k , in which Pk is the 

probability of each frequency being optimal in the brain wave (the total number is 512 and 

the probability is the optimal number of each frequency being divided by 512). Calculate 

the total entropy from the distribution of the probabilities of all 12 leads (total number is 

5 12* 12 )and the data is transmitted to the terminal processor E. 

[0154] 6 . Analysis module of single-frequency competition D6 

[0155] accumulate the number of the same optimal frequencies in frequency 

fluctuations f(n)of brain wave fluctuation signals along with time changing (1-512 

subsections ) so as to get the optimal probability curve. 

[0156] 7 . Analysis module of general pedigree of S pedigree D4a 

[0157] Form a graph of to present the general S pedigree signals produced by 

analysis module of S pedigree D4 and then the data and the graph are transmitted to the 

terminal processor E. 

[0158] 8 . Analysis module of fundamental pedigree D4b 

[0159] Analyze statistically the pedigrees of SK S2> S3, S4> S5> S6^ S7^ 
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S1K S13 which is corresponding to the nine fundamental frequencies of lmilli-Hz> 
2milli-Hz. 3milli-Hz> 4milli-Hz. 5milli-Hz. 6milli-Hz. 7milli-Hz> llmilli-Hz. 13milli-Hz 
in S pedigree respectively. When the frequency is higher than 3milli-Hz(3milli-Hz is 
included), the values of the multiple frequency should also be accumulated( for example, 
when 3milli-Hz is statistically analyzed, the values of 6mHz> 9mHz should also be 
accumulated,). Meanwhile, according to their positions in the brain, the 12 leads are 
divided into four sections to perform the statistical analysis, namely left anterior (F3 , F7, 
C3) , left posterior (T5, P3, 01). right anterior (F4, C4, F8)> right posterior (P4, 
T6 , 02 ) . The result is transmitted to the terminal processor E. 

[0160] 9. Analysis module of optimal value (D1-D8) D4c 
[0161] Present the power values of the dominant frequencies D1-D8 and their 
corresponding frequency in each lead according to the space position distribution of the 
lead and transmit the result to the terminal processor E. 

[0162] 10. Module for analysis of reversion of A/P and L/R unbalance D4d 
[0163] Said module D4d is to calculate the anterior-to-posterior A/P (such as 
F3/C3 , C3/P3 ) of the power value of each frequency according to the space distribution of 
the lead and then present the frequencies whose A/P values are more thanlO. Meanwhile, 
calculate the left-to-right ratio L/R (such as F3/F4, C3/C4) and present the frequencies 
whose L/R values are more than 10. 

[0164] 1 1 . Analysis module of characteristic spectrum line D4e 

[0165] The said module D4e includes analysis module of special frequency 
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D4ea 5 analysis module of different frequency D4eb, analysis module of continuum 
frequency D4ec and analysis module of optimal frequency D4ed, wherein: 

[0166] ( 1 ) analysis module of special frequency D4ea is to present the special 
frequency for each lead according to the space distribution position of the lead and the 
result is transmitted to the terminal processor E. 

[0167] (2) analysis module of different frequency D4eb is to present the 
different frequency for each lead according to the space distribution position of the lead 
and the result is transmitted to the terminal processor E. 

[0168] (3) analysis module of continuum frequency D4ec is to present the 
continuum frequency for each lead according to the space distribution position and the 
result is transmitted to the terminal processor E. 

[0169] (4) analysis module of optimal frequency D4edis to present the optimal 
frequency for each lead according to the space distribution position of the lead and the 
result is transmitted to the terminal processor E. 

[0170] 12. Analysis module of power space distribution of S pedigree D4f 

[0171] Select a certain spectrum line with a given frequency from the 
fluctuations of brain wave and present the power value of each lead according to the space 
position distribution of the lead. Then the power space distribution is transmitted to the 
terminal processor E. 

[0172] 13. Analysis module of power distribution D4g which includes analysis 
module of single- frequency power and distribution of the corresponding relative value 
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(L/R) D4ga, analysis module of the mean power distribution D4gb and analysis module 
of A/P, L/R D4gc. 

[0173] Add the power values corresponding to the optimal spectrum lines 
D1-D8 for each lead so as to get the total power value of each lead. The results are 
transmitted respectively to analysis module of single-frequency power and distribution of 
the corresponding relative value (L/R) D4ga, analysis module of the mean power 
distribution D4gb and analysis module of A/P, L/R D4gc. 

[0174] ( 1 ) analysis module of single-frequency power and distribution of the 
corresponding relative value (L/R) D4ga: present the values of fundamental frequency 
power and L/R ratio which is more than 10 or less than 0.1 according to the space 
distribution of the lead and the result is transmitted to the terminal processor E. 

[0175] (2) analysis module of the mean power distribution D4gb: present the 
mean power of each lead according to the space distribution of the lead and the result is 
transmitted to the terminal processor E. 

[0176] ( 3 ) analysis module of A/P, L/R D4gc: calculate the 
anterior-to-posterior ratio ( such as F3/C3 , C3/P3 ) and left-to-right ratio (such as F3/F4, 
C3/C4) of the power value according to space distribution of the lead and the result is 
transmitted to the terminal processor E. 

[0177] 14. Analysis module of long temporal dynamic curve D4h 

[0178] Said module D4h includes analysis module of the curve of long temporal 
S pedigree D4ha, analysis module of the curve of long temporal dominant spectrum lines 
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D4hb, analysis module of the curve of long temporal fundamental pedigree D4hc, analysis 
module of the curve of long temporal space distribution of power D4hd, analysis module 
of the curve of long temporal entropy D4he 5 analysis module of the curve of long temporal 
special frequency D4hf, analysis module of the curve of long temporal continuum 
frequency D4hg, analysis module of the curve of long temporal space distribution of 
fundamental pedigree D4hh, analysis module of the curve of long temporal conventional 
power spectrum D4hi 5 distinguishing module of long temporal event mark D4hj. 

[0179] Divide the long temporal (with sampling time longer than 18min) brain 
wave data into subsections according to the time length of 18 minutes and analyze the 
fluctuations of brain wave for each subsection (repeat the procedures from Dl to D4) . 
The result is transmitted to the following analysis modules. 

[0180] ( 1 ) analysis module of the curve of long temporal dominant spectrum 
lines D4hb is to form a curve by, using the frequency of the spectrum lines within the 
dominant spectrum line zone as the vertical axis and time as horizontal axis. A 
sub-interface or sub-window is provided in the display interface in the form of 
"sub-interface in interface" for arrangement selection of optimal lines (Dl-Dn) while 
displaying the curve of the curve of long temporal dominant spectrum lines. 

[0181] (2) analysis module of the curve of S pedigree D4hb is to form a curve 
by using fluctuation values of each pedigree or each spectrum line of each lead or all leads 
as vertical axis and time as horizontal axis so as to observe the changes of each frequency 
spectrum or pedigree in the whole brain along with time. A sub-interface or sub-window is 
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provided in the display interface in the form of "sub-interface in interface" for spectrum 
lines or pedigree selection while displaying the curve of S pedigree. 

[0182] (3) The analysis module of the curve of long temporal fundamental 
pedigree D4hc is to form a curve by using the fluctuation values in each time interval of 
the fundamental pedigree of all leads or each lead as vertical axis and time as horizontal 
axis so as to master the dynamic information of several fundamental pedigrees which have 
close relations to the brain's function. A sub-interface or sub-window is provided in the 
display interface in the form of "sub-interface in interface" for pedigree selection while 
displaying the curve of long temporal fundamental pedigree. 

[0183] (4)The analysis module of the curve of long temporal space distribution 
of power D4hd is to form a curve by using the power value of each spectrum line of 
each lead as vertical axis and time as horizontal axis so as to observe the changes of power 
value for each line in the lead along with time. A window is opened in the display 
interface in the form of "sub-interface in interface" for spectrum line selection, A 
sub-interface or sub-window is provided in the display interface in the form of 
"sub-interface in interface" for spectrum line selection while displaying the curve of long 
temporal space distribution of power. 

[0184] ( 5 ) analysis module on the curve of long temporal entropy D4he: Form 
a curve by using the entropy of each lead or all leads as vertical axis and time as 
horizontal axis, to present the changes of entropy along with time so as to observe energy 
distribution of the brain in different time. 
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[0185] (6) analysis module of the curve of long temporal special frequency 
D4hf: Form a curve by using the number of special frequencies of all leads as vertical axis 
and time as horizontal axis so as to observe the dynamic changes of special frequency 
along with time. 

[0186] (7) analysis module of the curve of long temporal continuum frequency 
D4hg: Form a curve by using the number of different frequencies in each lead or all leads 
as vertical axis and time as horizontal axis to observe decrease of patents' brain functions 
in different period. 

[0187] ( 8 ) The analysis module of the curve of space distribution of 
fundamental pedigree is to form curves each of which is presented according to the space 
position distribution of the lead by using power value of each lead which is acquired from 
the result of module for analysis of single-frequency power and the distribution of its 
corresponding relative value (L/R) D4ga as vertical axis and time as horizontal axis. A 
sub-interface or sub-window is provided in the display interface in the form of 
"sub-interface in interface" for pedigree selection while displaying the curve of space 
distribution of fundamental pedigree. So as to observe the changes of each fundamental 
pedigree along with time and space distribution. 

[0188] (9) the analysis module of the curve of conventional power spectrum 
D4hi is to form 8 dynamic curves by using power values of 8 frequencies D1-D8 which 
are arranged from big to small and selected from the conventional power spectrum which 
have maximal amplitudes and which are arranged in a descending order as vertical axis 
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and time as horizontal axis so as to observe the changes of optical frequencies in the brain 
wave along with time. 

[0189] (10) The analysis module of distinguishing of event mark D4hj is 
connected to each of the long temporal analysis modules respectively and it is applied to 
distinguish the event mark signals recorded in the brain wave recorder and present them in 
the playback of brain wave signals and the corresponding position on the time axis of 
various dynamic curves. 

Industrial Application 

[0190] According to methods and apparatus of the present invention, the data 
and curve, obtained by inspecting and calculating, can be used to analyze the situation of 
the brain function and neurotransmitter in the patents' brain. It can also be used to analyze 
the functional changes of the brain which is impossible for CT and NMR, and thus it can 
provide direct and objective index for the diagnosis of functional encephalopathy in the 
medical field and make up the blank of objective inspective index in functional 
encephalopathy such as lunacy in the medical field. 
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